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RAE ANTENNA-DAMPER PACKAGE 
DEPLOYMENT INTERACTION 


INTRODUCTION 

During the planning of the deployment sequence of the Radio Astronomy Ex- 
plorer, a gravity gradient stabilized spacecraft, a question arose as to when to 
deploy the libration damper package. The damper package is normally housed 
internal to the spacecraft and sometime during the spacecraft deployment phase , 
the package must be extended externally. The 24.5 pound package is accelerated 
by four one-pound negator springs to 2.3 ft/sec in traversing approximately 12 
inches from inside to outside the spacecraft and then is stopped abruptly. Two 
deployment times or methods were considered: the first was to deploy the pack- 
age prior to main boom deployment, during magnetic stabilization and the other 
was after main boom deployment. Thermal and dynamic effects tended to con- 
strain the first method, while the effect of damper package shock motion on the 
main booms constrained the other method. This analysis addresses itself to 
assessing the effect of deploying the damper package after main boom deployment. 


PROBLEM STATEMENT 

The orbital configuration of the RAE spacecraft with the "V" antennas de- 
ployed to an intermediate length of 450 ft is illustrated in Figure 1. 

Antenna deflections due to gravity and other orbital environment effects are 
neglected in this analysis. Since the deflections due to these effects are in the 
linear range, they can be added in a superposition fashion if desired. A free 
body diagram of a portion of one of the antennas is shown in Figure 2. The 
reaction acceleration, a 3 , experienced by the antenna boom during and after 
damper package deployment is resolved into components parallel and transverse 
to the antenna direction. This analysis considers the bending effect of the 
transverse loading but does not consider tension-compression effects which may 
result because of the parallel component of loading. The basic problem is to 
solve for the antenna deflection by adding the transverse loading terms due to the 
various accelerations to the standard linearized equation of beam bending: 


B 4 y 

3x 4 


acceleration loading per unit length 

El 


( 1 ) 


1 



4 ANTENNAS 
U2 in. dia., .002 i 
Be Cu 

p = 0.01404 Ib/ft 




Figure 2-Resolution of Reaction 
Acceleration into Parallel and 
Normal Components 

where 

x = distance along the axis of an undeformed antenna, ft (local coordinate 
system) 

y = deflection perpendicular to the undeformed antenna axis, ft 
El = flexural rigidity, lb -ft 2 
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r 


t 


The local coordinate system is embedded in the hub of the spacecraft and the 
hub is considered as a mass point. 

The acceleration loading consisted of three parts: 

1. That due to motion of antenna element with respect to local coordinate 
system: <' 2 v f >)t 2 

2. That due to lUoar motion of hub and hence local coordinate system as 
antennas deflect: d 2 ~ / ht 2 


3. That due to rigid body reaction motion of main spacecraft because of 
damper package deployment: a 3 

The first acceleration term,' 4 2 y '7n 2 , is self evident and needs no further 
explanation. The displacement of the hub, ^ , and hence ' Bt 2 can be 
expressed in terms of antenna deflection by using the principle of first mass 
moment. Referring to Figure 3, the first mass moment of the upper and lower 
antennas respectively is 



Upper V 



(2a) 



Figure 3- Antenna Coordinate System 
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where 


T? J - X cos $ 0 - y sin 0 Q (2b) 

and 

(2c) 
(2d) 


Lower V 
where 


f 

Vt\ 


V 2 ds 


t) zz x cos 0 n + y sin 


and s is a dummy variable of integration for x , Taking the difference and 
equating this to the product of total mass of spacecraft, M T , and center of mass, 
£ , we obtain 


or 


Wj. f = - 4 p sin 


£ = -4 P 


sin e 0 
M,. 



y(s, t) ds 


y(s, t) ds 


( 3 ) 


The acceleration is obtained by differentiation with respect to time twice 


3t 2 


4 p sin & Q 



it 2 


(s s t) ds 


( 4 ) 


Now the third acceleration term , a 3 , was determined experimentally together 
with Newton’s law of action and reaction as described in the next section of this 
report. 
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EXPERIMENTAL DETERMINATION OF a 3 

In obtaining a 3 as a function of time a deployment test of the damper pack- 
age was made under simulated zero-gravity conditions. The test gave the accel- 
eration experienced by the damper package. By an application of Newton 1 s third 
law, the reaction acceleration experienced by the rigid spacecraft hub and hence 
the booms is obtained. A schematic representation of the damper package as 
contained within the spacecraft is shown in Figure 4. There is a torsion wire 
supporting frame weighing approximately 3 lb which is driven downward in a 
track upon release until hard stops are reached. This occurs after a total travel 
of 12-3/8” at which time detents snap in place preventing back travel. Attached 
to this supporting frame through 4 shock absorbing springs (K - 70 lb/in.) is the 
main damper boom package weighing 21.5 lbs. This part of the damper package 
is free of its four teflon guide channels after about 4 inches of travel. 

The complete package assembly is driven by four negator springs of ap- 
proximately 1 lb force each. A maximum velocity of 2.3 ft/sec and a kinetic 
energy of 23.8 in. -lb are obtained by the time of impact with the hard stops. 

The four shock absorber springs absorb 11.4 in.- lb of the 21 it,.- lb contained 
in the main damper package by compressing .213” at which time they are bot- 
tomed. A rebound then begins which results in a reverse acceleration being 
applied. The force is applied to the detent pawls which are locked into place in 
the spacecraft structure. No shock absorption is provided in this direction so 
a hard rebound occurs here. All these phenomena are observed and recorded by 
an accelerometer mounted on the main damper package. 

An actual trace of the accelerometer output in g f s during a simulated damper 
package deployment is shown in Figure 5a. After an initial constant acceleration 
and then a deceleration pulse, the shock deceleration damps out completely in 
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approximately 1.5 seconds so only the initial and final portions of the curve are 
shown. In Figures 5b, 5c, and 5d, enlarged views of the initial portion of the 
curve and modifications performed upon it to get it into a form more suitable 
for analytical treatment are shown. In addition, brief explanations of the steps 
taken to simplify the curve are given. 

# 

The final form is represented analytically in two parts: the first part by a 
constant acceleration term and the second part by the product of an initial ac- 
celeration, an exponential decay function, and the sine of transient frequency 
times time. These parts are: 

r 

a = constant 0 < t < r (5a) 


a = a sin fit r ^ t <“ 

O 


(5b) 


From Newton's third law of action and reaction 



( 6 ) 


where m is the mass of the damper package and M is the mass of the spacecraft 
less damper package. Thus, when all of the acceleration loading terms are in- 
cluded in Equation 1, the governing differential equation for each antenna is of 
the form 

a4 y = __P_ 4p 2 sin 2 tf„ r 3 2 yfs,t)ds ^ a 3 sin d o n 

3x« El 3t 2 £IM t 3t 2 + El 

1 •'o 

where 

ma n Q 

a = 0 0 < t < r 

M t( A 2 + Q 2 ) 
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ACCELERATION IN g's ACCELERATION IN g's ACCELERATION IN g‘s 


1 - 


,9 SEC 


■^Yyyv 


"•* I 


- 1.5 SEC 


RELEASE 
OF PACKAGE 


(0) SHOCK WAVE EXPERIENCED BY 
MAIN DAMPER PACKAGE 



(b) ENLARGED VIEW OF INITIAL PORTION OF ACCELEROMETER 
READOUT 



A CONSTANT U MAY BE FOUND. 



Figure 5-Modificotion of the Transient Disturbance Wave for Analytical Treatment 
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and 


a = - - a n e’' Vv sin Q( t - r ) r <; t 

3 M 0 


The amplitude of the first pulse was adjusted so that the integral of the ac- 
celeration ( f ® a 3 d t ) is zero. This is a necessary condition so that no net 
velocity is imparted to the spacecraft due to the extension of the damper package. 

The boundary conditions are: 


and at t = 0 


y(0, t) = ^ (0, = o 

ox 

(8a) 

a^y(-C.t) _ djy(l.t) 

(8b) 

3x 2 3x 3 


y(x,0) =p (X ’ 0) =0 
at 

(8c) 


Equation 8a requires that the deflection and slope be zero at the origin. Equation 
8b requires that the moment and shear be zero at the end of the antenna and 
Equation 8c requires that the antenna be quiescent at t = 0. 

Equations 7 through 8c is an idealized dymanic statement of the antenna- 
damper package interaction. It is so formulated that only translation in the 
plane of the antennas is permitted to take place. Any rotational motion of the 
spacecraft hub due to misaligned deployment of the package is neglected. The 
governing differential equation will now be solved satisfying the specified boundary 
conditions and thus the deflection of a particular antenna and the resulting space- 
craft configuration change will be determined. 
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SOLUTION OF GOVERNING DIFFERENTIAL EQUATION 


The solution to Equation 7 can be iacilitated by expanding unity in terms of 
cantilever mode shapes X n (x) which satisfy the boundary conditions in Equations 
8a, 8b and are solutions of the fourth order differential equation: 


d< X n 

dx 4 



( 9 ) 


The X n (x) are dimensionless and normalized such that X ( i) = 1, and fi are the 
eigenvalues for the simple cantilever boundary conditions" (see references 1, 3). 

It is well known that the set of X (x) form an orthonormal vector basis such that 
an arbitrary function can be expanded in terms of these quantities. (See refer- 
ence 1.) Hence, 


00 


(X) 


( 10 ) 


where the constants C are given by the usual integration procedure 


A 

** A 


(x) dx 


r'i 

K 


± f'x, 

i J « n 


(x.) dx 


(x) dx 


Equation 11a can be evaluated by using Equation 9. Thus, 


(11a) 




X (x) dx = 

n N ' 


or 


% Jo " 


(x) dx = 


v w 

/ 3„ 4 


i 




X n (x) dx = - 


V (0) 


A 


(Hb) 


since there is no shear at the end of the antenna. Explicitly, the C is given by 

n 
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± V (0) 

n “ i ’<* 


( 12 ) 


n 


Equation 12 can be evaluated using tables given in reference 3. Substituting 
Equation 10 into Equation 7 results in the following 


B 4 


y + p ( 


ax 4 ei at 2 


B2y / 4p 2 sin 2 # C ^2 p a s l n 6 1 

7 ' 2 — 4 (s,t) ds + _ 3 

EIM,. Jo dt EI / i 


•)Z> 


(X) 


(13) 


Suppose we let the antenna deflection, y (x,t), have the following form 


00 


y(x, t) = V A n (t)X„(x) 


(14) 


where the A (t) which are functions only of t and have the dimension of length 
are to be determined. Substituting Equation 14 into 13 and rearranging results 
in 

' 4?..n-8 , r [' 




P 


Mr 


- J 


ds 


o 3t 2 

' C „ a 3 sin 6 0 I X „ ( x > = 0 


(15) 


Since each term in the parentheses is a function of time and the X n (x) are linearly 
independent, the coefficient of each X (x) must be zero. Thus, 


4p sin 2 # 
A + O) 2 A 2 

n n n 


Mj. 


I 




B 2 y (s, t) 
iT“2 


ds “ C n a 3 sin # 0 = 0 


n = 1, 2, 3 


(16) 
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where 


a? 


E^n 4 

P 


The integral in the above is evaluated by substituting Equation 14 into 16 and 
using Equation 11a, the definition of C R : 


»• 





pi 


s m‘ 


Qr 





a 3 sin d 0 


= 0 


(17) 


The infinite series in Equation 17 complicates the determination of the A , To 
circumvent this difficulty, a one term approximate solution will be obtained. 

That is, attention will be focused only on the determination of . Calculations 
not including the series term in Equation 17 showed that the first term contributed 
98% to the total deflection. With this approximation, the equation for A x is 


A, + a;? A, 


pi sin 2 0 o 

My 


CjAj - Cj a 3 sin d Q = 0 


or in more compact notation 


(18) 





+ 6J 2 Aj = Ba 3 


where 


(■- 


U)‘ 


p i sin 2 ^ 

mT 


o c 2 


) 


(19) 
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and 


B - 


sin 0 Q Cj 



P 


l sin 2 # 



) 


( 20 ) 


The denominator of Equations 19 and 20 indicates that pi sin 2 0 Q C 2 /M T must 
be much less than 1 for the approximation to be valid. Also, it will be noted that 
the natural frequency of the antenna system will be somewhat higher than the 
simple cantilever frequency. This is to be expected since the antenna system is 
tending toward a free-free beam frequency which is higher. For 0< t < r the 
solution for A x (t) satisfying the initial conditions: Aj ( 0) = Aj ( 0) = 0 (obtained 
from Equations 8c and 14) is 

B i 

A. (t) = — . (1 - cos co t) (21) 

co 2 


* 


where 



ma 0 QB 
Mr (H 2 + A 2 ) 


For r < t < oo the solution for A t ( t ) is 


B i B i 

Aj (t) = — sin ct)T sin (t — t) + — (1 - cos cor') cos a> (t - r) 


CO 


CO* 


( AG, - GG.\ 

-G 2 cos oj( t - r) + I-- — y s i n a> ( t -r) 


( 22 ) 


+ G 1 4 ~ T ) s in n ( t -r) +G 2 e"^( t “ T > cosfi(t - r ) 


where 


(A 2 + cu 2 -n 2 )B 2 

[(\ 2 + co 2 -n 2 ) 2 +4A 2 n 2 ] 


4 


(23) 
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( 24 ) 


2\HB 2 

' 2 [(\ 2 ' + « a - n 2 ) 2 + 4 \ 2 n 2 3 


where 


m a R 



It will be noticed that at t =r both solutions give the same value for \(r) and 
Aj (r) as they should. Further, it will be noticed that the exponential terms in 
Equation 22 rapidly damp out in a few seconds and so the deflection for the 
antenna several seconds after libration damper package deployment is essentially 
given by 


A t (t) 


B, 


or 


sin cor 




sin co (t - r) 


( 


B. 


+ | (1 - cos cot) - G, I COS CO (t -r) 


w 


(25) 


t > > r 


For wr small (which is true for RAE parameters). Equation 25 simplifies to 


V.» ■ ft 


B ,r X.G 0 fiG. . 

+ -- 1 sin co ( t - r) 

OJ 


B r 2 

+ l — !- G 2 ] cos co ( t - r) 


(26) 


Consistent with the one term approximation and the normalization discussed 
earlier, the tip deflection is determined by setting x = l in Equation 14 

/ y( 4 .t) = A 1 (t)X 1 ('i) = Aj(t) ( 27 ) 
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The hub movement is obtained from the following formula 


pisin & 0 Cj Aj (t) 

__ 


(29) 


which was derived from Equation 3. Another quantity of interest is the root 
moment acting on the antenna. This is given by the usual beam formula 


Moment = El = EIA (t) X! (0) 

Bx2 


(28a) 


= El 0 * Aj (t) 


(28b) 











DISCUSSION OF RESULTS 


Figure 6 shows a plot of antenna deflection versus time using the one term 
approximate solution. The numerical values of the physical parameters used 
are given in Appendix A. The salient points of the plot are, that maximum 
antenna deflection occurs 8 minutes after damper package deployment, and the 
maximum value is approximately 0.29 feet. The relatively long time to peak 
deflection is due to the inherently slow response of the long antenna element. 
The small deflection is due to the cancelling pulses of the damper package 
first accelerating and then decelerating. This can be seen by noting that for the 
RAE spacecraft, Cl is much larger than. k and w . With this approximation, 

G j andG 2 can be simplified to 


(k 2 + co 2 - Cl 2 ) B 2 


(30) 

[(\ 2 + co 2 -ci 2 ) 2 + 4A. 2 n 2 3 

n 2 + 3\ 2 

2 \nB 2 

2 \B 2 

(31) 

[(X 2 + co 2 - Cl 2 ) 2 + 4 A. 2 Cl 2 ] 

Cl(Cl 2 + 2k 2 ) 


where terms the order of k 2 /Ci 2 have been neglected compared to 1. Substi- 
tuting Equations 30 and 31 into Equation 26 and retaining terms the order of 
1/ co Cl results in the following expression for Aj ( t ): 


or 


Vt) 


a 0 mBQ 

M co 


1 


Cl 2 + k 2 


ci 2 + 3 k 2 


sin co(t- r) 


2a n mBH^ 2 sin co( t - r) 
Aj(t) = -2 — 

Ma; (Cl 2 + k 2 )(Cl 2 + 3 k 2 ) 


(32a) 

(32b) 


where 


B = 


sin 0 Q Cj 

' 1 -t£ sin2 


The ’’bucking” terms are seen in Equation 32 a and the net result is given in 
Equation 32b. Figure 7 shows the resulting changes of spacecraft configuration 
due to the deflections given by Equation 32b. Essentially the upper pair of 
antennas move outward while the lower pair move inward and vice versa in a 
sinusoidal fashion long after the applied transient load becomes small. This 
motion is one of the normal modes of vibration of the spacecraft identified in 
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reference 2. As was noted in the analysis, the antenna motion causes the hub of 
the spacecraft to move up and down about the equilibrium. Because of the 
small antenna deflections, the hub motion is small as well as the bending 
moment at the root of the antenna. 

Of concern during the RAE Program was the antenna deflection and whether 
the antenna motion damps out after a period of time. The vibration mode excited 
by the damper package deployment was uf'h! to initialize the RAE dynamics 
simulator and it was found that the antenna motion did not damp out but persisted 
indefinitely. This showed that the libration damper would not take out this type 
of antenna motion even if it was made functional a short time after deployment. 

There was apprehension as to how the RAE long antennas would behave in 
an orbital environment. A preliminary analysis of the damper package deploy- 
ment had indicated greater deflections of the long antennas than the results 
presented herein which were not complete at the time of finalizing the launch 
sequence. Because of the desirability of adhering to conservative approaches, 
it was decided to deploy the damper package prior to deploying the long 
antennas. 

In the actual flight of the RAE-A Spacecraft, the damper package was 
deployed successfully in this manner. 
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APPENDIX A 

NUMERICAL CALCULATIONS 
I. ANTENNA DEFLECTION 


Initial Values of Parameters 
= 20 n rad/sec 
a. - 3.37 sec 1 
a o ~ 6.54 g 
m = 24.5 lb 
M = 400 lb 
M t = 424.5 lb 
El = 2,130 lb in. 2 

p - 1.1697 x 10" 3 lb/in. or 3.03 x 10 “ 6 lb sec 2 /in. 2 
U Q = 30° 

g = 32.159 ft/sec 2 
i = 5400 in. 

in ^ 

<p x ( 0) - -1.46819 ^ 

> For Clamped-Free Beam (See Reference 3) 
= 1.8751041 J 

r = 0.9 sec 

Calculated Values 


1 . 




4 

1 
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2 2,13 x 10 3 / 1.8751041 V 

1 3.03 * 10’ 6 \ 5 4 * 10 3 / 


0.003196919 rad sec 


2. 


or 


4 K (Q) 

t , «4 



K (Q) 

2>ij 


upon substitution of X n - ; n 2 




1.46819 
^ 1 87 5104 1 ^ 


Cj = 1.565982 


(12)* 


in2 c i \ 
m t ) 

(0.003196919) 2 

/ (1.1697 x 10‘ 3 )(5.4 x 10 3 )(.5) 2 (1. 565982) 2 

\ ~ 424.5 

a 2 = 1.031439 x 10' 5 



(19) 


Numbers refer to equations in the main text 
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or 


* r 3.21602 x 10~ 3 


4. 


where 


2 a 


1 max 


Mu.(Q 2 + K 2 ) (fi 2 + 3-\ 2 ) 


B = 


sin d Q Cj 


t - 


p'l sin 2 t) Q C 2 


M. 


B = 


0.5 (1.565982) (32.159) 


1 - 


(1.1697 x 10- 3 ) (5.4 x 10 3 ) (.S) 2 (1.565982) 2 


B = 2S.41202 ft / sec 2 


A lmax = 2(6.54) ( 24. 5) ( 25 . 41 202 ) ( 2077 ) (3. 37) 2 

400(3.21602 x 10' 3 )(( 20 tt) 2 + ( 3. 37 ) 2 ) (( 20 77 ) 2 + 3(3. 37 ) 2 ) 

A , = 0. 286534 ft 

1 max 


A plot of deflection versus time is shown in Figure 6. The first maximum 
value of 0.287 ft occurs when t = 489.3 seconds. 
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